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ABSTRACT 

We have studied the rapid X-ray time variability in 99 pointed observations with the Rossi X-ray 
Timing Explorer (RXTE)'s Proportional Counter Array of the low-mass X-ray binary IE 1724-3045 
which includes, for the first time, observations of this source in its island and banana states, confirming 
the atoll nature of this source. We report the discovery of kilohcrtz quasi-periodic oscillations (kHz 
QPOs). Although we have 5 detections of the lower kHz QPO and one detection of the upper kHz QPO, 
in none of the observations we detect both QPOs simultaneously. By comparing the dependence of the 
rms amplitude with energy of kHz QPOs in different atoll sources, we conclude that this information 
cannot be use to unambiguously identify the kilohertz QPOs as was previously thought. We find that 
Terzan 2 in its different states shows timing behavior similar to that seen in other neutron-star low 
mass X-ray binaries (LMXBs). We studied the flux transitions observed between February 2004 and 
October 2005 and conclude that they are due to changes in the accretion rate. 

Subject headings: accretion, accretion disks — binaries: close — stars: individual (4U 1636-53,411 
1820-30,4U 1608-52,4U 0614+09,4U 1728-34, Terzan 2, IE 1724-3045) — stars: 
neutron — X-rays: stars 



1. INTRODUCTION 

Low-mass X-ray binaries (LMXBs) can be divided 
into systems containing a black hole candidate (BHC) 
and those containing a neutron star (NS). The accre- 
tion process onto these compact objects can be stud- 
ied through the timing properties of the associated X- 
ray emission (see, e.g.. | yan der K lis 2006, for a review). 
lHasinger fc van der Klisl (| 19891 ) classified the NS LMXBs 
based on the correlated variations of the X-ray spec- 
tral and rapid X-ray variability properties. They dis- 
tinguished two sub-types of NS LMXBs, the Z sources 
and the atoll sources, whose names were inspired by the 
shapes of the tracks that these sources trace out in an X- 
ray color-color diagram (CD) on time scales of hours to 
days. The Z sources are the most luminous, but the atoll 
sources are more numer ous and coyer a m uch wider range 
in luminosities (e.g. iFord et all 120001 . and references 
within). For each type of source, several spectral/timing 
states are identified which are thought to ari se from qual- 
itativ ely different inner flow configurations (jvan der Klisl 
20061). In the case of atoll sources, the three main states 
are the extreme island state (EIS), the island state (IS) 
and the banana branch, the latter subdivided into lower- 
left banana (LLB), lower banana (LB) and upper banana 
(UB) states. The EIS and the IS occupy the spectrally 
harder parts of the color color diagram and correspond 
to lower levels of X-ray luminosity (L x ). The associated 
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patterns in the CD are traced out in hours to weeks. 
The hardest and lowest L x state is t he EIS, which shows 
strong (up to 50% rms amplitude, see lLinares et al.|[2007l 
and references within) low-frequency flat-topped noise 
also known as band- limited noise (BLN) . The IS is spec- 
trally softer and has higher X-ray luminosity than the 
EIS. Its power spectra are characterized by broad fea- 
tures and a dominant BLN component which becomes 
weaker and generally higher in characteristic frequency 
as the flux increases and the > 6 keV spectrum gets 
softer. In order of increasing L x we then encounter the 
LLB, where twin kHz QPOs are generally first observed, 
the LB, where 10-Hz BLN is still dominant and finally, 
the UB, where the < 1 Hz (power law) very low fre- 
quency noise (VLFN) dominates. In the banana states, 
some of the broad features observed in the EIS and the 
IS become narrower (peaked) and occur at higher fre- 
quency. In particular, the twin kHz QPO can be found 
in the LLB at frequencies higher than 700 Hz (the lower 
of them with frequencies generally lower th an a 1000 Hz, 
while the upper one up t o 1258 Hz, see e.g.. lBelloni et al.l 
[20051 : 1 Jonker eTa!]l2007h . only one kHz QPO can be gen- 
erally found in the LB, and neither of them is de t ected in 
the UB (see reviews by Ivan der Khsll2000l 12004 I2006D . 

A small number of weak NS LMXBs which do not get 
brighter than a few times 10 36 ergs s _1 (usually burst 
sources and often referred t o as 'weak' or 'faint bursters', 
see e.g, iMuno et al.l 120031 and references within) re- 
semble atoll sources in the EIS, but in the absence of 
state transitions this identification has been tentativ e 
(see for example iBarret et al J 120055 Ivan der Klisll2006l) . 
An important clue is provided by the correlations be- 
tween _Jie_conmonejitfr 

e.g. Ivan Straaten et all l2002t 120031 : lAltamirano et all 
2008) which helps to identify components across sources. 
For example, Ivan Straaten etafl ([2001 [2003) compared 
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the timing properties of the atoll sou rces 4U 0614+09, 
4U 1 608-52 and 4U 1728-34 (see also lAltamirano et all 
l2008l for similar results when the atoll source 4U 1636- 
53 was included in the sample) and conclude that 
the frequencies of the variability components in these 
sources follow the same pattern of correlations when plot- 
ted versus the frequency of the upper kHz QPO {v u ). 
Ivan Straaten et all (|2003t ) also showed that low luminos- 
ity systems extend the frequency correlations observed 
for the atoll sources. This last result gave further clues in 
t he link between the atoll and th e low luminosity sources. 

iPsaltis fc Chakrabartvi (|1999f ) found an approximate 
frequency correlation involving a low-frequency QPO, 
the lower kHz QPO frequency and two broad noise com- 
ponents interpreted as low-frequency versions of these 
features. This correlation spans nearly three decades 
in frequency where the Z and bright atoll sources pop- 
ulate the > 100 Hz range and black holes and weak 
NS systems the < 10 Hz ran ge. As already noted by 
IPsaltis fe Chakrabartvi (|1999l ). because the correlation 
combines features from different sources which show ei- 
ther peaked or broad components with relatively little 
overlap, the data are suggestive but not conclusive with 
respect to the existence of a single correlatio n covering 
this wide frequency range (|van der Klisll200l . 

The low-luminosity neutron star systems can play a 
crucial role in clearing up this issue. Observations of 
different source states in such a system could connect 
the < 10 and > 100 Hz regions mentioned above by 
direct observation of a transition in a single source. 
In the case of the patter n of correlations reported by 
Ivan Straaten et ahl (|2003f ). low luminosity NS systems 
extend the frequency correlations observed for ordinary 
atoll sources down to ~ 100 Hz. Unfortunately, the low 
luminosity NS systems are usually observed in only one 
state (EIS), which makes it difficult to properly link these 
sources to the atoll sources. However, some of these 
objects show rare excursions to higher luminosity lev- 
els which might correspond to other states. The oc- 
currence of these excursions are usually unpredictable. 
Therefore, in practice it was not possible until now to 
check on the frequency behavior of the different variabil- 
ity components as such a source enters higher luminosity 
states. 

IE 1724-3045 is a classic low luminosity LMXB; a per- 
sistent Low-Mas s X-ray binary locate d in the globular 
cluster Terzan 2 (|Grindlav et al.lfl980f ) which is a metal- 
rich globular cluster of the galactic bulge. Its distance is 
estim ated to be between 5.2 to 7.7 kpc (jOrtolani et al.1 
Il997fl . These values are consistent with that derived from 
a ty pe I X-ray burst that s howed photosp heric expansion 
(see Grindlav et al .11 19801 but also see iKuulkers et al.l 
12003. iGallowav et alj 1 20061 ) . The type I X-ray bursts 
observed from this source also indicate that the compact 
objec t is a weakly magnetize d neutron star dSwank et al.l 
[19771 iGrindlav et alTfl980l) . lEmelvanov et all (120021 ) 



have shown, using ~ 30 years of data from several X- 
ray satellites, that the luminosity of Terzan 2 increased 
until reaching a peak in 1997, after which it started to de- 
crease. They suggest that the evolution of the donor star 
or the infl uence of a th i rd sta r coul d be the caus e of this 
behavior. lOlive et alj (|1998l ) and iBarret et all (|2000al) 
have shown that during earlier observations of Terzan 2 
its X-ray variability at frequencies > 0.1 Hz resembled 



that of black hole candidates. This state was tentatively 
identified as the extreme island state for atoll sources. 
Until now, no kilohertz quasi-periodic oscillations have 
been reported for this source, which was attributed to 
the fact that t he source was always observed in a single 
intensity state (jBarret et al.ll2000H ). 

Monitoring observations by the All Sky Monitor 
aboard the Rossi X-ray Timing Explorer showed that the 
source was weakly variable in X-rays (less than about a 
factor of 3 on a few day time scale for the first 8 years of 
the m onitoring). However. recentlv lMarkwardt fc Swankl 
(2004) reported (u sing PCA monitoring obser vations of 
the galactic bulge - ISwank fc Markwardtfe OOl) that dur- 
ing February 2004, IE 1724-3045 flared up from its rel- 
atively steady ~ 20 mCrab to ~ 66 mCrab (2-10 keV). 
In this paper we report a complete study of the timing 
variability of the source. For simplicity, and since only 
one bright X-ray source is detected in the globular clus- 
ter (See Section 14. 1|) . in the rest of this report we will 
refer to IE 1724-3045 as Terzan 2. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Light curves and color diagrams 

We use data from the Rossi X-ray Timing Explorer 
(RXTE) Proportional Counter Array (PCA; for instru- 
ment information see IZhang et al.l 119931 ; iJahoda et al.l 
2006). There were 534 slew observations until October 
30 t?l 2006, which are part of the PCA monitoring obser- 
vations of the galactic bulge (|Swank fc Markwardtj|200"l[) 
and which were performed typically every 3 days. These 
observations were only used to study the long-term L x 
behavior of the source. 

There were also 99 pointed observations in the nine 
data sets we used (10090-01, 20170-05, 30057-03, 50060- 
05, 60034-02, 80105-10, 80138-06, 90058-06 & 91050-07), 
containing ~ 0.8 to ~ 26 ksec of useful data per observa- 
tion. We use the 16-s time-resolution Standard 2 mode 
data to calculate X-ray colors. Hard and soft color are 
defined as the 9.7-16.0 keV / 6.0-9.7 keV and 3.5-6.0 keV 
/ 2.0-3.5 keV count rate ratio, respectively, and intensity 
as the 2.0-16.0 keV count rate. The energy-channel con- 
version was done using the pca_e2c_e05v02 table provided 
by the RXTE Team 5 . Channels were linearly interpo- 
lated to approximate these precise energy limits. X-ray 
type I bursts were removed, background was subtracted 
and deadtime corrections were made. In order to correct 
for the gain changes as well as the differences in effective 
area between the PCUs themselves, we normalized our 
colo rs by the correspond i ng Crab Nebula color va lues; 
(see IKuulkers et al.lll994t Ivan Straaten et alj 120031 see 
table 2 in lAltamirano et al Eoosf for average colors of the 
Crab Nebula per PCU) that are closest in time but in the 
same RXTE gain ep och, i.e., with the sa me high voltage 
setting of the PCUs (jJahoda et al.ll2006D . 

2.2. Fourier timing analysis and fitting models. 

For the Fourier timing analysis we used either the 
Event modes E_125us_64M_0.1s or E_16us.64M_0.8s, or 
the Good Xenon data. Leahy-normalized power spec- 
tra were constructed using data segments of 128 seconds 
and 1/8192 s time bins such that the lowest available 

5 see http://heasarc.gsfc.nasa.gov/docs/xte/pca_news.html 
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frequency is 1/128 » 8 x 10 -3 Hz and the Nyquist fre- 
quency 4096 Hz. No background or deadtime correc- 
tions were performed prior to the calculation of the power 
spectra. We first averaged the power spectra per obser- 
vation. We inspected the shape of the average power 
spectra at high frequency (> 2000 Hz) for unusual fea- 
tures in addition to the usual Poisson noise. None were 
found. We then subtracted a Poisson noise spectrum 
estimated from the power bet ween 3000 and 40 00 Hz, 
using the method developed bvlKlein-Woltl (12004D based 
on the analytical function of I Zhang et al.l ( 19951 ). In this 
frequency range, neither intrinsic noise nor QPOs are ex- 
pected based on what we observe in other sources. The 
resulting p ower spectra were converted to squared frac- 
tional rms (|van der Klis 1995). In this normalization the 
power at each Fourier frequency is an estimate of power 
density such that the square root of the integrated power 
density equals the fractional rms amplitude of the intrin- 
sic variability in the source count rate in the frequency 
range integrated over. 

In order to study the behavior of the low-frequency 
components usually found in the power spectra of neu- 
tron star LMXBs, we needed to improve the statis- 
tics. We therefore averaged observations which were 
close in time an d had both similar color s and power 
spectra (see e.g. Ivan Straa ten et al.l 1200 2. 2003, 2005; 
lAltamirano et all I2005L l2008|. and referen ces within - 
see also Appendix in lAltamirano et al.ll2008l for a discus- 
sion on other possible methods). 

To fit the power spectra, we used a multi-Lorentzian 
function: the sum of several Lorentzian components 
plus, if necessary, a power law to fit the very low fre- 
quency noise (VLFN - see Ivan der K lis 2006] for a re- 
view). Each Lorentzian component is denoted as Lj, 
where i determines the type of component. The char- 
acteristic frequency {v max as defined below) of Li is 
denoted z/j. For example, L u identifies the upper kHz 
QPO and v u its characteristic frequency. By analogy, 
other components go by names such as Lg (lower kHz), 
L[ w (lower Lorentzian) , L^Hz (hectohertz), (hump), 
Lb (break frequency), L^ (second break frequency) and 
their frequencies as Vi, V£ ow , VhHz, Vh and Vb, v b2, re- 
spectively. Using this multi-Lorentzian function makes 
it straightforward to directly compare the characteris- 
tics of the different components observed in Terzan 2 to 
those in previous works which used the same fit function 
(e.g..lBelloni et alJl20ljljyan Straaten et alJ 12001 120051 : 
lAltamirano et aill2005l 120081 . and references therein). 

Unless stated explicitly, we only include those 
Lorentzians in the fits whose single trial significance ex- 
ceeds 3<7 based on the error in the power integrated from 
to oo. We give the frequency of the Lorentzians in 
terms of characteristic frequency v max as introduced by 
iBelloni et all (|200l : v max = + (FWHM/2) 2 = 

+ 1/(4<3 2 ). For the quality factor Q we use the 
standard definition Q = is /FWHM. FWHM is the full 
width at half maximum and the centroid frequency of 
the Lorentzian. The quoted errors use A% 2 = 1.0. The 
upper limits quoted in this paper correspond to a 95% 
confidence level (A% 2 = 2.7). 



2.3. Energy spectra 



Since the energy spectra of the quiet state (see Sec- 
tion 13. ip of Terzan 2 have alrea dy be en studied in previ- 
ous wo rks (see e.g. lOlive et~aLlll998l: iBarret et al.lll999l 
l2000al) . in this paper we concentrate on the 14 obser- 
vations that sample the flaring period (see Section l3~T|) . 
In all 14 cases, we used data of both the PCA and the 
HEXTE instruments. 

For the PCA, we only used the Standard 2 data of PCU 
2, which was active in all observations. The background 
was estimated using the PCABACKEST version 6.0 (see 
FTOOLS). We calculated the PCU 2 response matrix for 
each observation using the FTOOLS routine PCARSP 
V10.1. For the HEXTE instrument, spectra were accu- 
mulated for each cluster separately. Dead time correc- 
tions of both source and background spectra were per- 
formed using HXTDEAD V6.0. The response matrices 
were created using HXTRSP V3.1. For both PCA and 
HEXTE, we filtered out data recorded during, and up 
to 30 minutes after passage through the South Atlantic 
Anomaly (SAA). We only use data when the pointing 
offset from the source was less than 0.02 degrees and the 
elevation of the source respect to the Earth was greater 
than 10 degrees. We did not perform any energy selec- 
tion prior to the extraction of the spectra. Finally, we 
fitted the energy spectra using XSPEC V11.3.2L 




50000 50500 51000 51500 52000 52500 53000 53500 54000 
Time (MJD) 

Fig. 1. — Above: PCA count rate obtained from the monitoring 
observation s of Terzan 2 (Monit oring observation of the galactic 
bulge - see Swank & Markwardt 2001). These data was used for 
the study of the long-term variability (see Sections 12. 4l fc l3^8t . Be- 
low: Terzan 2 intensity (Crab normalized - see Section O versus 
time of all pointed observations. These observations were use d for 
the st udy of the 0.1—1200 Hz X-ray variability (see Sections 12.21 
[3741 & EO]} except for the observation 90058-06-04-00, which is 
marked with an arrow (see Section [3). The modified Julian date is 
defined as MJD = Julian Date-2400000.5. 



2.4. Search for quasi-periodic variations over days and 

months 

Recently, IWen et ail (|2006f ) have performed a system- 
atic search for periodicities in the light curves of 458 
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Time interval 
(MJD) 


Duration 
\ (lavs ) 


Maximum Count Rate 
(cts/s/5PCU) 


Flare 
label 


N°. of pointed obs. 
sampling the flare 


~ 53038 - 53071 


33 


~ 850 


Fl 


7 


~ 53127 - 53155 


28 


~ 950 


F2 


4 


~ 53233 - 53250 


17 


~ 830 


F3 


2 


~ 53493 - 53502 


9 


~ 450 


F4 





~ 53566 - 53602 


36 


~ 1150 


F5 


1 


~ 53631 - 53651 


20 


~ 650 


F6 





~ 53934 - 53972 


38 


~ 865 


F7 






TABLE 1 

Data on the 6 flares observed until MJD 53667. See Section for details. (The modified Julian date is defined as MJD = 

Julian Date-2400000.5). 



sources using data from the RXTE All Sky Monitor 
(ASM). Terzan 2 was not included in their analysis, prob- 
ably due to the fact that the ASM source average count 
rate is low: 2.05 ±0.01 count/s (the average of the errors 
- l/nJ2 err i ~ is 0.8 count/s). 

Since the PCA gal actic bulge monitoring 
(| Swank fc Markwardtl [2001) has observed the source 
for more than 8 years, and the lowest detected source 
count rate was 170 ± 5 counts/sec, this new data set 
provides useful information to search fo r long term 
modulations. Lomb-Scargle pe riodograms (|Lomblll97a 
IScargle|[l982t iPress et aHll992D as well as t he phase dis- 
persi on minimization technique (PDM - see Stellingwcrl 
were used. The Lomb-Scargle technique is ideally 
suited to look for sinusoidal signals in unevenly sampled 
data. The phase dispersion minimization technique is 
well suited to the case of non-sinusoidal time variation 
covered by irregularly spaced observations. 

3. RESULTS 
3.1. The light curve 

Figure [1] shows both the PCA monitoring lightcurve 
of th e source (see upper panel and lSwank &: Markwardtl 
I2001D and the Crab normalized intensity (see lower panel 
and Section 12. ip of each pointed observation versus 
time (in units of modified Julian date MJD = Julian 
Date— 2400000.5). In the rest of this paper, we will refer 
as "quiet period" to that between MJD 51214 and 52945, 
and as "flaring period" between MJD 52945 and 53666. 

During the quiet period, 333 monitoring measurements 
of the source intensity and 85 pointed observations sam- 
ple the behavior of the source. The count rate slowly de- 
creases from an average of ~ 300, to an average of ~ 190 
counts/s/5PCU at an average rate of —0.059 ± 0.002 
count/s/day. In the flaring period, 7 flares sampled 
with 201 monitoring observations were detected with the 
galactic bulge scan. 14 pointed observations partially 
sampled parts of 4 of these flares. In Table [1] we list ap- 
proximate dates at which the flux transitions occurred, 
the flare durations and the maximum count rates de- 
tected with the PCA. As mentioned in Section 12. 4i the 
monitoring is done approximately once every three days; 
additional gaps in the data are present due to visibility 
windows. As of course we do not have details of flares 
that may have occurred during these gaps, the informa- 
tion in Table[T]is only approximate. In Figure[2]we show 
the intensity of the source during the flaring period. We 
label the different flares Fl, F2, F3, F4, F5, F6 and F7 
in order of time of occurrence. 

We detected three Type I X-ray bursts. One was dur- 
ing the quiet-state observation 10090-01-01-021 and two 



during the flaring-state observations 80138-06-06-00 and 
90058-06-02-00. A detailed study of these X-ray bursts as 
well as a compa rison to bursts observe d in other sources 
can be found in lGallowav et all (|2006f ). 




53000 53200 53400 53600 53800 54000 
Time (MJD) 



Fig. 2. — Above: Intensity in units of the Crab Nebula of the 
observations which were performed during the flare state. The 
different states are labeled (see Section 13.21 1. Arrows mark the 
observatio ns in which kHz QPOs were detected (see also Section l3.3l 
and Table l3~3)l . Below: Intensity versus time showing the seven 
flares. The data are the same as those of Figure [T] 

3.2. Color diagrams; identification of states 

Figure [3] (top) shows the color-color diagram of the 99 
pointed observations. For comparison, we also include 
the color-color diagram of the atoll source 4U 1608-52, 
which has been observed in all extreme island, island 
and banana states (Figure [3l bottom) . The similarity in 
shape suggests that Terzan 2 underwent state transitions 
during observations of the flares. Based on Figure [3] 
we can identify the probable extreme island, island 
and banana state with the hardest, intermediate and 
softest colors, respectively. Since partially sampled 
patterns in the color color di agrams are not nec essarily 
unambiguous, (see review by Ivan der Klid 12006). power 
spectral analysis (below) is required to confirm these 
identifications. The extreme island state is sampled 
with 85 pointed observations which are clumped in 
2 regions at similar hard colors but at significantly 
different soft colors. We find 3 observations in the 
island state. They sampled the lowest luminosity 
sections of flares F2, F3 and F5 (see Figure [2]). The 
banana state is sampled by 11 pointed observations: 
7 during Fl, 3 during F2 and 1 during F3. The 
identifications above are strengthened by the similar- 
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ities in power spectral shapes bet ween Terzan 2 
those reported 



and 



2003 
2002 
200E 



other sources ( van Straaten et al 



Di Salvo ct al 



Di Salvo et al 



iLinares et alJ 




van Straaten ct al 



van Straaten et al 
20051 : 

Migliari et all 120051 : lAltamirano et all l2008f ). In the 
following sections we describe the power spectra in more 
detail. 




1.1 1.15 1.2 1.25 1.3 1.35 1.4 
Soft color (Crab) 




EIS 



M 1 1.1 1.3 1.3 1.4 

Soft color (Crab) 

Fig. 3. — Top: Hard color versus soft color normalized to Crab as 
explained in Section [2] Different symbols represent the selections 
used for averaging the power spectra as explained in Section [2] and 
shown in Figure [6] (See Figure [T] for symbols). The arrow marks 
observation 90058-06-04-00, which was excluded from interval N 
(see Section [3}. Bottom: Hard color versus soft color normal- 
ized to Crab for the NS source 4U 1608-52. This source has been 
observed in all expected atoll states: extreme island state (EIS), 
island state (IS), lower left banana (LLB), lower banana (LB) and 
upper banana (UP). The similarity between both figures suggest 
that Terzan 2 has been observed in similar states as 4U 1608-52. 



3.3. kHz QPOs 

We searched each averaged observation's power spec- 
trum for the presence of significant kHz QP Os at frequen- 
cies x *, 400 Hz. As r eporte d in other works (|Barret et al.l 
119991 iBelloni et all 12001 ). during the quiet period the 
power spectra of single averaged observations show sig- 
nificant power up to ~ 300 Hz which is fitted with 
broad (Q < 0.5) Lorentzians plus if necessary, one sharp 
Lorentzian to account for L^p. During the flaring pe- 
riod, we found that several observations show power ex- 
cess above 400 Hz. For each observation we fitted the 



averaged power spectra between 400 and 2000 Hz with 
a model consisting of one Lorentzian and a constant to 
take into account the QPO and Poisson noise, respec- 
tively. In 5 of the 14 observations that sample the flaring 
period, we detect significant QPOs, with single trial sig- 
nificances up to 6.5cr. In Table [3~3l we present the results 
of our fits and information on these observations. As can 
be seen, the first three observations were performed dur- 
ing the rise of the first flare while the fourth and fifth 
observations where done during the decay of the second 
flare (see Figure [2]). As can be seen in Table [3~3l there 
are significant frequency variations. Unfortunately, due 
to the sparse coverage of the flares and the fact that we 
cannot detect the QPO on shorter time scales than an 
observation, no further conclusions are possible. 

We do not significantly detect two simultaneous kHz 
QPOs in any of the 5 observations. In order to search for 
a possible second kHz QPO, we used the shift -and-add 
method as described bv lMendez et al.l (|1998bf j. We first 
tried to trace the detected kilohertz QPO us ing a dynam- 
ical p ower spectrum (e.g. see figure 2 in iBerger et al.1 
1996) to visualize the time evolution of the QPO fre- 
quency, but the signal was too weak to be detected on 
timescales shorter than the averaged observation. There- 
fore, for each observation we used the fitted averaged 
frequency (see Table 13. 3|) to shift each kilohertz QPO 
to the arbitrary frequency of 770 Hz. Next, the shifted, 
aligned, power spectra were averaged. The average power 
spectrum was finally fitted in the range 300-2048 Hz so 
as to exclude the edges, which are distorted due to the 
shifting method. To fit the averaged power spectrum, we 
used a function consisting of a Lorentzian and a constant 
to fit the QPO and the Poisson noise, respectively. We 
studied the residuals of the fit, but no significant power 
excess was present apart from the 770 Hz feature. In 
Figure |4] we show the fitted kHz QPO for observation 
80138-06-02-00 (no shift and add was applied) and the 
shifted-and-added kHz QPO detected with the method 
mentioned above. Since it is known that t he kHz QPOs 
become stronger at higher energies (e.g. IBerger et al.1 
119961: iMendez et al"1ll998at Ivan der Klisll2000L and refer- 
ences within), we repeated the analysis described above 
(which was performed on the full PC A energy range), us- 
ing only data at energies higher than ~ 6 keV or higher 
than ~ 10 keV. Again, no significant second QPO was 
present. It is important to note that this method can 
produce ambiguous results as we cannot be sure we are 
always shifting the same component (either L u or Lg). 
We also tried different subgroups, i.e. adding only two 
to three different observations, but found the same re- 
sults. 

To investigate the energy dependence of the kHz 
QPOs, we divided each power spectrum into 3, 4 or 5 en- 
ergy intervals in order to have approximately the same 
count rate in all the intervals. We then produced the 
power spectrum as described in Section [2] and refitted 
the data where both frequency and Q were fixed to the 
values obtained for the full energy range (see Table l3~3|) . 
In Figure[5]we show the results for the representative kHz 
QPOs in flares Fl and F2 (observations 80138-06-03-00 
and 90058-06-02-00, respectively). Similarly to what is 
observed in other sources, the fractional rms amplitude 
of the kHz QPOs increases with energy. The data show 
that there is no significant difference in the energy de- 
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MJD 


Obsld 


Used in 
interval 


vq (Hz) 


FWHM 


Q 


rms (%) 


Sig. 


c/s 
(5PCU) 


53054.94 


80138-06-02-00 


M 


731 ±3 


33.2 ± 7.3 


22.1 ±4.8 


10.4 ±0.8 


6.5cr 


~ 665 


53055.00 


80138-06-02-01 


M 


764 ± 2 


14.3 ±5.2 


53.4 ± 19.4 


6.9 ±0.8 


4.3cr 


~ 690 


53056.69 


80138-06-03-00 


M 


772 ± 1 


15.1 ±4.2 


51.1 ± 14.2 


7.1 ±0.6 


5.6cr 


~ 660 


53145.57 


90058-06-01-00 


N 


559 ±2 


17.3 ±6.9 


32.3 ±12.8 


6.6 ±0.7 


4.7cr 


~ 765 


53147.21 


90058-06-02-00 


N 


599 ± 1 


20.3 ±4.8 


29.5 ±6.9 


7.4 ±0.5 


6.6cr 


~ 610 



TABLE 2 

Fit results for the 5 averaged observations showing significant kHz QPOs. 



pendence of the kHz QPO at ~ 599 Hz with that at 
~ 772 Hz. 




500 1000 

Frequency (Hz) 



Fig. 4.— Fit to the kHz QPO in observation 80138-06-02-00. In 
the subplot on the top-right of the figure, we show the fit to the 
data after using the shift-and-add method on all 5 observations 
reported in Table [3751 where the ma in pe ak was set to the arbitrary 
frequency u = 770 Hz (see Section 13-36 - No significant detections 
of a second kHz QPO were found. 



3.4. Averaged power spectrum 

As mentioned in Section [2~2l in order to study the be- 
havior of the low-frequency components usually found in 
the power spectra of neutron star LMXBs, we needed 
to improve the statistics. We therefore averaged obser- 
vations which were close in time and had both similar 
colors and power spectra. The resulting data selections 
are labeled from A to R (ordered mainly in time - see Ta- 
ble [3] for details on which observations were used for each 
interval and their colors). Their corresponding average 
power spectra are displayed in Figure [Sj 

3.4.1. The power spectra of the quiet period 

Intervals A to L are all part of the quiet segment. 4 to 
6 components were needed to fit all the power spectra of 
this group, where 11 out of the 12 power spectra showed 
significant broad components at ~ 150, ~ 10, ~ 1 and at 
~ 0.2 Hz. Interval G is the exception where the broad 
component at ~ 150 Hz was not significantly detected 
(8.5% rms-amplitude upper limit). A QPO (Q > 2) at 
~ 0.8 Hz was also significantly detected in 11 out of 




4U 0614+09 A 
■ ■ [ i h 4U 1608-52- Bergeretal. ^-o-^ 

1 10 
Energy (keV) 

Fig. 5. — Energy dependence of two representative kHz QPOs 
of Terzan 2. One in flare Fl (pentagons - OBSid 80138-06-03- 
00) and another in flare F2 (downward pointed triangles - OBSid 
90058-06-02-00). For comparison we show the energy dependence 
of both lower (circles & open squares ) and upper ( squares) peak 
of the atoll source 4U 1608-52 berger et al. 1996; Mcnd ezeTaT] 
1998b) and the upp er peak (upwards po i nted triangles) of the ato ll 
source 4U 0614+09 l IMendez et al.H1997l ; lvan Straaten et aI.II2000P ). 

12 power spectra, interval B being the exception (1.4% 
rms-amplitude upper limit). Finally, an extra component 
(L v i , where vl stands for "very low" ) at frequency v v i ~ 
0.1 Hz was detected only in interval C. 

Similar power spectra have been reported in the ex- 
treme islan d state of the atoll sou rces 4U 0614+09, 
4U 1728-34 (Ivan Straaten et al.ll2 002). XTE J1118+480, 
SLX 1735-269 (IBelloni et all I2002T) and 4 U 16 08-52 
(Ivan Straaten et all l2003f) . lOlive et all (|1998l ) and 
iBelloni et al.l (|2002h analyzed an early subset of the data 
we present in this work. Our results are consistent 
wit h the frequencies r e ported in those works. F ollow- 
ing Ivan Straaten et al.lj2002t ). IBelloni et all (|2002f ) and 
Ivan Straaten et ail ([20031 ) ■ we identify the components 
as L u , L tow , L h , L LF and L b , where v u ~ 150 Hz, 
viow ~ 10 Hz, Vh ~ 1 Hz, Vlf ~ 0.8 Hz and vi> ~ 0.1 Hz, 
respectively. These identifications are strengthened by 
the correlations shown in Figure [5J in which we plot the 
characteristic frequency of all components versus that of 
v u for several atoll sources (see Section l3.6.1[) ; clearly, we 
observe the same components as seen in the EIS in other 
atoll sources. In Figure [7] we provide a modeled power 
spectra were the main 5 Lorentzians (i.e. L u , Li OWl Lh, 
Llf and Lb) are represented. Components in Figures [6] 
and [7] are represented with the same line style. 

3.4.2. The power spectra of the flaring period 

Intervals M, N, O, P, Q and R are part of the flaring 
period. For intervals M, N and Q, high frequency (> 
400 Hz) single QPOs are seen which can be identified as 
either the lower or the upper kHz QPO. 
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Frequency (Hz) 

Fig. 6. — Power spectra and fit functions in the power spectral density times frequency representation. Each plot corresponds to a different 
region in the color-color and color-intensity diagrams (see Figures [3] and [T} . The curves mark the individual Lorcntzian components of the 
fit. The re-binning varies between the panels to balance the frequency resolution with the signal-to-noise ratio. For a detailed identification, 
see Table H 



Interval A 



Observation 



Soft color (Crab) 



Hard color (Crab) 



Intensity (Crab) 



10090-01-01-000 
10090-01-01-001 
10090-01-01-00 
10090-01-01-020 
10090-01-01-021 
10090-01-01-022 
10090-01-01-02 



1.3615 ± 
1.3508 ± 
1.3619 ± 
1.3783 ± 
1.4435 ± 
1.3638 ± 
1.3723 ± 



0.0017 
0.0019 
0.0019 
0.0021 
0.0019 
0.0017 
0.0039 



1.0760 
1.0761 
1.0757 
1.0781 
0.9846 
1.0788 
1.0801 



±0.0013 
±0.0014 
±0.0015 
±0.0016 
±0.0013 
±0.0013 
± 0.0029 



0.0389 ± 
0.0384 ± 
0.0384 ± 
0.0386 ± 
0.0468 ± 
0.0387 ± 
0.0389 ± 



0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 



Interval B 



20170- 


05- 


01 


00 


1 


3501 


± 





0075 


1 


0660 


± 





0058 





0383 


± 





0001 


20170- 


05- 


02 


00 


1 


3440 


± 





0077 


1 


0492 


± 





0059 





0384 


± 





0001 


20170- 


05- 


03 


00 


1 


3498 


± 





0076 


1 


0541 


± 





0058 





0389 


± 





0001 


20170- 


05- 


01 


00 


1 


3449 


± 





0073 


1 


0632 


± 





0056 





0394 


± 





0001 


20170- 


05- 


05 


00 


1 


3492 


± 





0071 


1 


0703 


± 





0055 





0398 


± 





0001 


20170- 


05- 


06 


00 


1 


3392 


± 





0072 


1 


0650 


± 





0056 





0394 


± 





0001 


20170- 


05- 


07 


00 


1 


3465 


± 





0072 


1 


0726 


± 





0055 





0393 


± 





0001 


20170- 


05- 


08 


00 


1 


3295 


± 





0071 


1 


0605 


± 





0056 





0382 


± 





0001 


20170- 


05- 


09 


00 


1 


3398 


± 





0070 


1 


0704 


± 





0054 





0389 


± 





0001 


20170- 


05- 


10 


00 


1 


3163 


± 





0064 


1 


0687 


± 





0051 





0387 


± 





0001 


20170- 


05- 


11 


00 


1 


3601 


± 





0071 


1 


0714 


± 





0054 





0400 


± 





0001 


20170- 


05- 


12 


00 


1 


3365 


± 





0074 


1 


0747 


± 





0058 





0399 


± 





0001 


20170- 


05- 


13 


00 


1 


3532 


± 





0070 


1 


0721 


± 





0053 





0409 


± 





0001 


20170- 


05- 


11 


00 


1 


3215 


± 





0076 


1 


0768 


± 





0061 





0407 


± 





0001 


20170- 


05- 


15 


00 


1 


3417 


± 





0070 


1 


0658 


± 





0054 





0394 


± 





0001 


20170- 


05- 


16 


00 


1 


3398 


± 





0081 


1 


0669 


± 





0062 





0398 


± 





0001 


20170- 


05- 


17 


00 


1 


3290 


± 





0072 


1 


0832 


± 





0057 





0386 


± 





0001 


20170- 


05- 


18 


00 


1 


3318 


± 





0086 


1 


0682 


± 





0067 





0392 


± 





0001 


20170- 


05- 


19 


00 


1 


3336 


± 





0068 


1 


0667 


± 





0053 





0395 


± 





0001 


20170- 


05- 


20 


00 


1 


3309 


± 





0069 


1 


0599 


± 





0054 





0397 


± 





0001 



TABLE 3 

Observations used for the timing analysis. The colors and intensity are corrected by dead time and normalized to the 
Crab Nebula (see Section [2j- The complete table can be obtained digitally from Ap,I. 
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0.01 0.1 1 10 100 1000 10 4 



Frequency (Hz) 

Fig. 7. — Example of the typical power spectra of the quiet 
period. Th e main 5 Lorentzians are L u , Li ow , L^, L^p and L\> 
(see Section 15, 4. 

The kHz QPOs in interval M and N correspond to the 
averages of significant QPOs observed in single observa- 
tions. As seen in Table 13.31 the characteristic frequency 
of the kHz QPOs averaged in each of the two intervals are 
within a range of 50 Hz. By our averaging method we are 
affecting the Q value of the kHz QPOs but improving the 
statistics for measuring the characteris tics of the features 
at low er frequencies (see Appendix in lAltamirano et all 
2008, for a discussion on this issue). 

Interval Q is an average of three single observations 
(see Table [3]) that individually do not show significant 
QPOs ( Q > 2) at high frequencies although low-Q power 
excess can be measured. in interval Q is only 2.6a 
significant (single trial) but required for a stable fit. 

Interval O shows a broad component at 30.5 ± 6.3 Hz 
and a 3.3% rms low-frequency noise. In this case, 
the very low frequency noise was fitted with a broad 
Lorentzian because a power law gave an unstable fit. 
The excess of power at v > 1000 Hz is less than 3a 
significant. Interval P shows a power spectrum with a 
power-law low frequency noise, and two Lorentzian com- 
ponents at frequencies 15.5 ± 1.6 and 30.5 ± 6.3 Hz (see 
Figure [5]). In this case the high x 2 /dof = 218/163 re- 
veals that the Lorentzians do not satisfactorily fit the 
data. As can be seen in Figure O there is a steep de- 
cay of the power above v ~ 35 Hz and power excess at 
~ 70 Hz. A fit with three Lorentzians becomes unsta- 
ble. To further investigate this, we refitted the power 
spectrum using instead two Gaussians and one power 
law to fit the power at v < 40 Hz, and one Lorentzian 
to fit the possible extra component at ~ 70 Hz. The 
steeper Gaussian function better fits the steep power de- 
cay than the Lorentzians. In this fit, with three more 
free parameters, we obtain a x 2 /dof — 188/160, and the 
Lorentzian at 69.1 ±2.6 Hz becomes 3.4c single-trial sig- 
nificant. This power spectrum i s very simi lar to those re- 
ported bv lMigliari et alj (|2004 12005) and iMigliari et all 
(2005) for the atoll sources 4U 1820-30 and Ser X-l, re- 
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Fig. 8. — The characteristic frequencies u max of the var- 
ious power spectral components plotted versus u u . The grey 
symbols mark the ato ll sources 4U 0614+09, 4U 1728-3 4 
(van Straatc rTetafll2003). 4U 1608-52 llvan Straaten et alj|2003h, 
Aql X-l |Reig et al.l "l2004l and 4U 1636-53 lAltamirano et al.l 
2008) and the low luminosity bu rsters Terzan 2 (previous results), 
GS 1826-24 and SLX 1 735-269 Ivan Straaten et al.ll2005l, but also 
see IBell oni et al. 2002). The black bullets mark our results for 
Terzan 2. Note that we only plot the result s for Intervals A-N and 
Q (note that as mentioned in Section l3.6.1l the points for intervals 
M and N (y u > 800 Hz) were plotted under the assumption that 
v u = ft + 300 Hz). For i ntervals O, P and R no kHz QPOs were 
detected (see Section l3,4jl . 

spectively. Besides the similarity in shape of the power 
spectrum, it is interesting to note that in both cases these 
authors found a best fit with components at similar fre- 
quencies to the ones we observe in Terzan 2 and which 
they interpreted as Vb, Vh an d v hHz- This coincidence 
suggests that the sources were in very similar states. To 
our knowledge, no systematic study of this state has been 
reported as yet. 

Interval R consists of only one observation (90058-06- 
05-00) of ~ 1.4 ksec of data. In addition to a power law 
with index a = 3.1 ± 0.7, we detect one Lorentzian at 
68.5 ± 7.2 Hz. Its frequency is rather high if we compare 
it with the other power spectra presented in this paper 
and even when compared with results in other sources 
(see Figure EJ. This result might be due to blending 
of components due to the low statistics present in this 
power spectrum. 

From the 99 pointed observations, only observation 
90058-06-04-00 was not included in any of the averages 
described above. The averaged colors of this observation 
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are similar to those of Interval N (observations 90058- 
06-01-00, 90058-06-02-00) but the power spectrum does 
not show a significant QPO at ~ 560 Hz and can be 
fitted with a single Lorentzian at v max — 23.6 ± 3.3 Hz, 
Q = 0.5 + 0.2 and rms= 13.2 ±0.9%. The residuals of the 
fit show excess power at ~ 800 Hz but no significant kHz 
QPO. Although the colors of observation 90058-06-04-00 
are similar from those of interval N, the power spectrum 
is different. Because of this, we refrained from averaging 
90058-06-04-00 into interval N. 




0.85 0.9 0.95 
Hard Color (Crab) 

Fig. 9. — Fractional rms (%) amplitude versus hard color normal- 
ized to the Crab Nebula as explained in Section [2] The light-grey 
squares represent the data of the "quiet" state while the dark-grey 
circles represent the data for the "flare" state. The black circles 
represent the upper limit (Ax 2 = 2.7) for the observations 80138- 
06-05-00 and 80138-06-06-00 which also partially sampled the flare 
states. 



plot the characteristic frequency of all components ver- 
sus that of v u for the atoll sources 4U 0614+09, 4U 1728- 
34, 4U 1608-52, 4U 1636-53 and Aql X-l and the two 
low l uminosity bursters GS 1826-2 4 and SLX 1735- 
269 (Ivan Straaten et ap 12001 120031 : iReig et al.l 120041 : 
lAltamirano et all [2008). Using these correlations to 
identify the highest frequencies we observe (in inter- 
vals M, N & Q) as either v u or V(, presents a prob- 
lem. Both interpretations give consistent results since 
the correlations observed are complex when v u > 600 Hz 
(|van Straaten et al.ll2005l ). 



3* * 



PBK99-BPK02 
Terzan 2 (This paper) i 



Fig. 10. — PBK relatio n after Psaltis & Chakrabarty (PBK99 
- [l9§9Tl and IBelloni et all fBPK02~ T2002l) . The black squares at 
vlf < 10 Hz represent the data of Intervals A and C-L while the 
black squares at vlf > 20 Hz represent the data of Intervals M 
and N. In Interval B we do not detect Llf and in intervals Q-R 
we do not detect Lf. 



3.5. Integrated power 

In order to study the rms amplitude dependence on 
color and intensity, we calculated the average integral 
power per observation between 0.1 and 1000 Hz. In Fig- 
ure[H]we show the 0.01-1000 Hz averaged rms amplitude 
(%), of each of the 99 observations, versus its average 
hard color. The observations which sample the island 
and banana state correlate with the averaged hard color. 
This type of correlation has been already observed in 
black hole candida tes and neutron star systems (see e.g. 
iHoman et al.l l2001). At colors harder than 1.0, there are 
two clumps (grey squares in Figure [H]) which can also be 
seen in the color-color diagram (Figure [3|) . They both 
correspond to observations of the extreme island state 
(quiet state) of the source. As we show in Figure O the 
power spectral shape remains approximately the same 
with time. However, for a given hard color, the total rms 
amplitude can change up to ~ 30%. No correlation with 
time, intensity or soft color was found and these changes 
are seen within a small range in intensity (less than 4 
mCrab). Similar results are also observed when the rms 
amplitude is calculated in the 0.01-300 Hz range. 

3.6. Comparing Terzan 2 with other LMXBs 

3.6.1. The flaring period 

In FigureQJJwe plot the frequency correlations b e tween 
Llf and Lg rep orted b y Psalti sfc Chakrabartvi (1999) 
and updated by IBelloni et al.l (j2002fl . In Figure [8] we 



In recent work. IBarret et ahl ([2005a, be) have system- 
atically studied the variation of the frequency, rms am- 
plitude and the quality factor Q of the lower and upper 
kHz QPOs in the low-mass X-ray binaries 4U 1636-536, 
4U 1608-52 and 4U 1735-44. Although Q depends on the 
frequency of the component, these authors show that Qi 
is always above ~ 30, while Q u is generally below ~ 25. 
By comparing our re sults for the kHz QPOs with those of 
IBarret et all (|2005allbl lcl). we find that the 5 QPOs listed 
in Table |3"U1 (and averaged in intervals M and N) are rel- 
atively high-Q and hence consistent with being Lg, but 
also still consistent with being L u . The quality factor of 
the kHz QPO found in interval Q is low (2.4 ± 0.6) and 
hence this kHz QPO is probably L u . 

In Figure [TJJ we plot the fractional rms amplitude 
of all components (except Llf) versus v u for the 4 
atoll sources 4U 0614+09, 4U 1728-34, 4U 1608-52, 
4U 1636-53. With respect to the kHz QPO identifi- 
cation, the interpretation that we found Lg in interval 
N is strengthened by the rms amplitudes of the two 
low-frequency components found in the averaged power 
spectrum. If the QPO we observe is not Lg but L u , 
then the low frequency QPO pairs can be identified as 
either L^-Lf, or L^-L^i- For v u — 586.1 ± 6.6, the 
pair Lh-Lf, is not consistent with what we observe for 
other atoll sources (Figure [TT]), since Lb is not seen with 
rms amplitude as low as 3.9^ g%. The pair L^-L^ 
is not consistent with the data either, since L{,2 is al- 
ways observed at v u > 800 Hz. If the QPO is Li, then 
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based on what we obse r ve in other well-studied s ources 
(IMendez et all I1998U IMendez fc van der Klisl 119991: 
Pi Salvo et"aLl 120031: IBarret et all I2005U Iva n der Klisl 
20061 IMendez fc Bellonil l2007f) we expect that the fre- 
quency difference between kHz QPOs is A„ = v u — ~ 
300 and therefore v u ~ v t + 300 ~ 586 + 300 ~ 886 Hz. 
Under this assumption (i.e. ~ 300) and by using 
the same reasoning as above, then only the pair L^-L^ 
is consistent with the data. In the case of interval M, 
only one component is found at low frequencies with 
an rms amplitude of 6.4 ± 0.6%. This result is consis- 
tent with several interpretations when compared with 
the data shown in Figure [TT1 Therefore, for interval M 
we cannot improve confidence in the identification of the 
kHz QPO using Figure HU 

In Figures [5] [TO] & HH we have plotted the data for 
Terzan 2 based on the identifications above. As discussed 
later in this paper, such identifications need to be con- 
firmed. 

3.6.2. The quiet period 

In Figure [Tl] we show that the rms amplitude of 
L u , Lg OWl Lh and in Terzan 2 approximately 
follow the trend observed for other sources. Since 
iPsaltis fc Chakrabartvl (|1999t ) interprets Li ow and Lg as 
the same component in different source states, in Fig- 
ure QT] we plot the data for both components together. 
Of course, in Figure [11] and Figure [8J as well, there 
is the well-known gap between these two components. 
Regarding lower frequency components, the point in- 
side the circle represents our result for L v i in interval 
C, which is the weakest component found in the EIS 
of Terzan 2. We plotted our results with those for 
Lb2- This com ponent might be related with the VLFN 
Lorentzian (see | Schnerr eta l. 2003; Rccri nk et al . 1 120051) . 
lAltamirano et al. (2008) showed that the rms amplitude 
of Llf of several atoll sources does not correlate with v u . 
There is also no evidence for a correlation in the case of 
Terzan 2 (not plotted). 

3.7. Luminosity estimates from spectral fitting 

We fitted the PCA and the two cluster's HEXTE spec- 
tra simultaneously using a model consisting of a black 
body to account for the soft component of the spec- 
tra and a power law to account for the hard compo- 
nent. In some cases, it was necessary to add a Gaus- 
sian to_take_jjito accou nt the iron Ka line (6.4 keV, see 
e.g. iWhite et~all[l98l . We ignore energies below 2.5 
and above 25keV for the PCA spectra, and below 20 
and above 200 keV for the HEXTE spectra (see e.g. 
IBarret et al.l l2000a| ). In most cases, it is not possi- 
ble to well constrain the interstellar absorption nu if 
we lack spectral information below 2.5 keV. We there- 
fore opted to fix nn to the value nn — 1.2 ■ 10 22 
H atoms cm " 2 (in the Wisconsin cross s ection wabs 
model - see iMorrison fc McCammonl 119831 ) based on 
previous ASCA/BeppoS AX results (|01ive etal] 119981 : 
IBarret et al]ll999L l2000ah . 

Assuming a distance of 6.6 kpc, we found that all 
14 observations have luminosities between ~ 0.4 and 
~ 1.35 x 10 37 erg s _1 in the energy range 2-20 keV. 
We also found that at high energies (20-200 keV) the lu- 
minosities of most observations were less than 0.09 x 10 37 
erg s _1 . The exceptions are the three observations which 



sample the island state, which show 20-200 keV lumi- 
nosities of ~ 0.16, ~ 0.23 and ~ 0.29 x 10 37 erg s _1 
(observations 90058-06-03-00, 90058-06-06-00 and 91050- 
07-01-00, respectively). This may be compared with ob- 
servations of the brightest interval of the quiet period of 
Terzan 2, which have averaged luminosities Li_20fceV = 
0.81 x 10 37 erg s" 1 a nd L 2 o-2oofc e y = 0.48 x 10 37 erg s" 1 
(jBarret et al.ll2000af) . Clearly, in between the flares the 
luminosity can drop to similarly low values as in the quiet 
period. This is consistent with the lightcurve we show 
in Figure [H We no te that the observations studied by 
IBarret et al.l (|2000aD correspond to MJDs 50391 - 50395 
(November 4-8, 1996) and sample the brightest part of 
the quiet period of this source observed with RXTE (see 
Figure [J) . 

We are particularly interested in the luminosity at 
which the kHz QP Os are detec t ed in Terzan 2 compared 
with other sources. iFord et ail (]2000f ) have measured si- 
multaneously the properties of the energy spectra and 
the frequencies of the kHz QPOs in 15 low-mass X-ray 
binaries covering a wide range of X-ray luminosities. The 
observations of intervals M and N (see Table [T]) have av- 
erage luminosities L2-50keV /L e dd between 0.025 and 0.04 
(where L e dd — 2.5 x 10 38 erg s -1 ). The three observa- 
tions that sample the island state (interval Q) have aver- 
age luminosity L2-50fceV / Ledd — 0.02. This means that 
during the flares, Terzan 2 shows kHz QPOs at similar 
luminosities to the atoll sources Aql X-l, 4 U 1608-52, 
4U 1 702-42 and 4U 1728-34 (see figure 1 in IFord et all 
[2003). 

3.8. Lomb Scargle Periodograms 

During the quiet period (51214-52945 MJD), we found 
no significant periodicities using either the Lomb-Scargle 
or the PDM techniques in the full data set nor in sub- 
intervals. 

As shown in Figure [U the flares seem to occur ev- 
ery ~ 60 - 100 days. Both Lomb-Scargle and the PDM 
techniques confirm this with a significant signal of pe- 
riod P ~ 90.55 days. In Figure [12] we show the PCA 
lightcurve (top) versus a 20-bin 90.55 days period folded 
lightcurve (bottom). The folded lightcurve matches the 
occurrence of most of the flares. However, it is clear 
that the flares are not strictly periodic. For example, F3 
seems to occur later and F6 occur earlier than expected. 
Furthermore, it is not possible to say if F7 is an early or 
late flare, or even a blend of two flares (we observed a 
small flare which peaked at ~ 53926 MJD, followed by a 
big one which peaked at ~ 53958 MJD). Although there 
are gaps in the data, Figure [TO] suggests that some flares 
do not occur at all (see arrow in this figure). 

4. DISCUSSION 

4.1. Contamination by a second source in the same field 

of view? 

As shown in Figure [TJ the luminosity of the source 
slowly decreases with time during the quiet period 
51214-52945 MJD. Although the rms amplitude changes 
up to 30% (see Figure \§§ , the X-ray timing characteris- 
tics are very similar (see Interval A to L in Figure [6]). 
During the 53000-53700 MJD period, the source shows 
flares which show different X-ray timing characteristics 
consistent with the island and banana states observed 
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Fig. 11. — The fractional rms amplitude of all components (exce pt L^p) plotted versus i/ u - The symbols are labeled in the plot. The 
data for 4U 1728-34, 4U 1608-52 and 4U 0614+09 were taken from [van Straat en~et~aTI (20051 ). The data for 4U 1636-53 were taken from 
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L v i while the points inside the square represent results for L; om (see Section [4] for a discussion). Note that as mentioned in Section 13.6.1 1 
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Fig. 12. — The PCA monitoring observation lightcurve of 
Terzan 2 in the time interval 53013-54050 MJD (above) versus 
a series of folded light curves with period of 90.55 days (below). 
The bl ack a rrow shows the position of a possible missed flare (see 
Section Ol 



in other atoll sources (see e.g. Ivan Straaten et al.1 [20031 
120051 [Belioni et al.ll200l lAltamirano et all 1200.1 120081) . 
A possible mechanism of the observed flux variations in 
Terzan 2 could be the emergence of a second X-ray source 
in this globular cluster unresolved by the 1° (FWHM) 
field of view of the PCA. If two sources are observed si- 
multaneously with RXTE, then we would expect to see 
power spectra which are a combination of the intrinsic 
time variability of both sources. 

To further investigate this, we compared the absolute 
rms amplitude that we observe both in the quiet and 
flaring states. Observation 80105-10-01-00 is the last ob- 
servation performed during the quiet period from which 
we measured an average source countrate of ~ 200 counts 
/ second. The integrated power between 7.8 • 10~ 3 and 
1 Hz is (3.4±2)10 -2 , which corresponds to a fractional 
rms amplitude of ^ 18 ± 0.6%, i.e. an absolute rms am- 
plitude of 36 ± 1 counts/second. 

Observation 80138-06-01-00 is the first observation per- 
formed during the flaring state. Its average countrate is 
~ 740 counts / second and the absolute rms amplitude 
in the 0.0078 — 1 Hz range is 21.4 ± 4.4 counts/second 
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(2.9 ± 0.6% fractional rms amplitude). Clearly, the 
absolute rms amplitudes are different when comparing 
quiet and flaring periods. If we repeat the analysis 
using the second RXTE observation during the flaring 
state (80138-06-02-00), the discrepancy is higher. This 
observation has an average source countrate of ~ 405 
counts/second and the upper limit for the absolute rms 
amplitude in the 7.810~ 3 — 1 Hz frequency range is 5 
counts/second. 

Given the characteristics of the power spectra, flares 
cannot be explained by assuming that another source 
has emerged, unless Terzan 2 turned off at the same time 
that the other X-ray source turned on, which is unlikely. 
Therefore, we conclude that the flux transitions are in- 
trinsic to the only low mass X-ray bina ry detected in the 
globu lar cluster Terzan 2: IE 1724-30 (|Revnivtsev et all 
12001 . 

4.2. The kilohertz QPOs, different states and their 
transitions 

The results presented in this paper show that the low 
luminosity source IE 1724-3045 in the globular cluster 
Terzan 2 can be identified as an atoll source. This is the 
first time a sou rce previously class i fied as weak burst 
sourc e (see e.g, iBelloni et al.l l2002t Ivan Straaten et all 
l2003l and references within) showed other states than 
those of the extreme island state, confirming previous 
suggestions that these sources are atoll sources. We have 
identified the new states as the island and banana states 
based on comparisons between color color diagrams of 
different sources and the characteristics of the power 
spectra. We have detected at least one of the the kHz 
QPOs, and as explained in Sections 13.4.21 and 13.31 in 5 
cases we may be detecting the lower kHz QPO (intervals 
M and N - see also Table I3.3[) and in one case the upper 
one (interval Q). No simultaneous twin kHz QPOs were 
detected within any of the 14 observations that sample 
the flares. Future observations of flares will allow us to 
confirm these identifications and might allow us to detect 
both kHz QPOs simultaneously. 

We found that the frequencies of the various compo- 
nents in the power spectra of Terzan 2 followed previ- 
ously reported relations (Figures [8] and [10]) . Terzan 2 is 
a particularly important source in the context of these 
frequency correlations because it is one of the few neu- 
tron star sources that has been demonstrated to show 
power spectral features that reach frequencies as low 
as ~ 0.1 Hz, which is uncommon for neutron star low 
mass X-ray binaries, but not for black holes. Our re- 
sults demonstrate that in each of the flares, Terzan 2 un- 
dergoes flux transitions that, if directly observed, would 
probably allow us to resolve current ambiguities in the 
identification of components, such as the case of Lg ow 
component in atoll sources. This component is inter- 
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quencies (see 
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2002) which becomes peaked 



at higher frequencies, while other authors interpret Lg 
and Le ow as different compo nents (see e.g. discussion 
in Ivan Straaten et ail l2003f ) . Another example is the 
ident ification of the upper kH z QPO at low frequen- 
cies. Ivan Straaten et alX(|2003l ) have suggested that the 
broad component observed at ~ 150 Hz in the EIS 
of atoll sources becomes the peaked upper kHz QPO 



L u . These authors based their interpretation on the 
fr equency correlations show n in Figure [5] Nevertheless, 
as Ivan Straaten et alj (|2003l) argue, these identifications 
should be taken as tentative. One way to confirm the link 
between them would be to observe the gradual transfor- 
mation from one to another one. 

During the time between flares, the source shows inten- 
sities similar to those measured before the quasi-periodic 
flares started. Unfortunately there are no observations 
during those intervals, but we expect that then Terzan 2 
shows X-ray variability similar to that reported in inter- 
vals A-L. If this is the case, the state transition between 
the extreme island state and the island state should be 
observable in observations at the beginning or at the 
end of each flare. Given the relatively gradual and pre- 
dictable transitions, Terzan 2 becomes the best source 
known up to now to study these important transitions. 

4.3. On the ~ 90 days flare recurrence 

The quasi-periodic variations over days and months 
observed in some LMXBs X-ray light-curves are gener- 
ally associated with the possible precession period of a 
tilted accretion disk or alternatively long term periodic 
variations in the accretion rate or periodic outbursts of 
X-ray transients. Some examples are the ~ 35 cycle in 
Her X-l which is thought to be caused by a varying ob- 
scuration of the neutron star by a tilted-twisted process- 
ing accretion disc; the ~ 1 70 days accretion cycle of the 
atoll source 4U 1820-30, (jPriedhorskv fc Terrell H984at 
ISimonl 120031 ): the 122-125 day cy cle in the outbursts of 
the recurrent transient Aql X-l (jPriedhorskv fe Terrelll 
ll984bUKitamoto et al.lll993l ). Understanding the mech- 
anisms that trigger the long-term variability associated 
with variations in the accretion rate of LMXBs can allow 
us to better predict, within each source, when the state 
transitions occur. This is useful because these transitions 
are usually fast and therefore difficult to observe. 

The power spectra of our observations of Terzan 2 dur- 
ing the flaring confirm that the source undergoes EIS-IS- 
LLB-LB-UB state transitions, as observed in other neu- 
tron star atoll systems (and not as see n for Z-sources, 
see reviews by Ivan der Klis 2004 l2006l and references 
within). As the source increased in X-ray luminosity, 
we found that the components in the power spectra in- 
creased in frequency which is consistent with the inter- 
pretation that the accretion disk is moving inwards to- 
ward the compact object. Therefore, the flaring with av- 
erage 90 days period is most probably an accretion cycle. 
We note that the modulation of the light-curve could be 
related to the orbital period of the system or set by the 
precession of a tilted disk. However, the mechanisms in- 
volved in those interpretations are very unlikely to affect 
the frequency of the kHz QPOs. 

If the flares are explained as an accretion cycle, then 
it is puzzling why the source underwent a smooth de- 
crease of L x for ~ 8 years before it started to show the 
flares. Terzan 2 may not be the only source that shows 
this kind of behavior. For example, KS 1731-260 is a 
low-luminosit y burster that has show n a h igh L x phase, 
during which iRevnivtsev fe Sunvaevl (|2003l ) reported a 
possible ~ 38 days period, and a low L x phase, during 
which much stronger variability was observed (which was 
described as red noise). After i ts low L x phase, KS 173 1- 
260 has turned into quiescence (jWiinands et alj2002bllah . 
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In Figure [13] we show the bulge scan light curve of the 
source during the low L x phase. At MJD ~ 51550 the 
source reached very low intensities, then flared up again 
for > 250 days to finally turn into quiescence. The low 
luminosities are co nfirmed by t he AS M light curve (not 
plotted). Recently. IShih et all ((2005f t reported that the 
persistent atoll source 4U 1636-53 has also shown a pe- 
riod of high L x followed by a period of low L x . During 
high L x , no long-term periodicity was found, but a highly 
significant ~ 46 days period was observed after its L x de- 
cline. 

These similar patterns of behavior might point towards 
a common mechanism, which then must be unaffected by 
the intrinsic differences between these sources. 

For example, while Terzan 2 remained with ap- 
proximately constant luminosity in its extreme island 
state for ~ 8 years before showing long term period- 
icities, 4U 1636-53 and KS 1731-260 were observed 
with variable luminosity and in different states, includ- 
ing the banana state in which the k Hz QPOs were 
found (see e.g. IWiinands et al.1 Il997l and IShih et al. 
Bool for 4U 1636-53 andlWiinands k. van der Klid ll997 
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and iRevnivtsev fc Sunvaevi 120031 for KS 1731-260) 
While Terzan 2 reached a maximum luminosity of 
L x /LEdd — 0.02 during one of the flares, 4U 1636- 
53 shows similar luminosities only at its lowest L x 
levels (while it has r eached L x /LEdd ^ 0.15 - see 
lAltamirano et all [2008). Further differences may be re- 
lated to whether these systems are normal or ultra- 
compact binaries. While 4U 163 6 -53 i s not ultra- 
compact (see below) . lin't Zand et al.l (|2007f ) has recently 
proposed that Terzan 2 may be classified as ultra- 
compact based on measurements of its persistent flux, 
long burst recurrence times and the hard X-ray spectra. 
If the luminosity behavior of these sources is related, the 
differences outlined above suggest that the mechanism 
that triggers the modulation of the light curve at low L x 
may not depend on the accretion history the luminosity 
of the source or even whether the system is ultra-compact 
or not. The modulation period may depend on these fac- 
tors. 

Unfortunately, we cannot compare the orbital periods 
and the companions of the three systems, as these are 
only known for 4U 1636-53 (~ 3.8 hours and ~ 0.4 M 
ICasares et al.|[20 06). Nevertheless, with the present data 
it is already possible to exclude some mechanisms. For 
example , mass trans fer feedback induced by X-ray irra- 
diation fpsaki 1985)) is unlikely. In this model, X-ray 
radiation from the compact object heats the companion 
star surface, causing enhancement of the m ass acc retion 
rate in a runaway instability. However, in [Qsakif s sce- 
nario, it is not clear how the system could remember 
the phase of the cycle if one of the flares is missed or if 
the size of the flares differs much. Flares F4 and F5 in 
Terzan 2, independently of the other two sources, may 
already raise an objection to this model. Although we 
miss part of F4 due to a gap in the data, Figure Q] shows 
that F4 was quite short (less than 9 days), while F5 was 
t he longest (5, 36 d ays) and strongest flare. 

Shih et all (|2005f ) have suggested that the atoll source 
4U 1636-53 may turn into quiescence after its low L x 
period, as was observed for KS 1731-260. Such an ob- 
servation for 4U 1636-53 as well as for Terzan 2 would 
give credibility to the link between these sources. To our 
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Fig. 13. — The PCA monitoring observation lightcurve of the 
atoll source KS 1731-260 during part of its low L x period. Unfor- 
tunately, there are no PCA monitoring observations of the source 
before to MJD 51200. Clearly, the source flares up similarly to 
Terzan 2 before it turns into quiescence. Interestingly, the data 
at MJD ~ 51550 shows that the source had a period of very low 
intensity, followed by a flaring up that lasted for ~ 250 days before 
the source finally turned into quiescence. 

knowledge, there is no model which predicts such behav- 
ior. 

4.4. Energy dependence as a tool for kHz QPO 
identification 

iHoman fc van der Klis! ((2000) discovered a single 
695 Hz QPO in the low mass X-ray binary EXO 0748-676 
and identified this QPO as the lower kHz QPO. These 
authors based their identification on the fact that at that 
time: (i) from the 11 kilohertz QPO pairs found in atoll 
sources, eight had ranges of lower peak frequencies that 
include 695 Hz, which was the case for only three of the 
upper peaks and (ii) the upper peaks in atoll sources 
generally had Q lower than ~ 14, while their QPO had 
Q> 38, value more common for lower peaks. While from 
Figure [5] it can be seen that (i) is not valid anymore, 
since the upper kHz QPOs have been detected down to 
300 -400 H z (and possibly down to < 100 Hz - see Sec- 
tion 13.4. ip , at these low frequencies L u is usually much 
broader than they observed, w hich confirms their identi- 
fi cation (see Sectionl3.6.1l andlBarret et al.l l2005al lbl[cl) . 

Ho man fc van der Klisl (|2000l ) also based their iden- 
tification on the comparison of the energy dependence 
of the QPO with that of the two kilohertz peaks in 
4U 1608-52, which have rather different energy de- 
pendencies (the power-law rms-ene rgy relation f or Li 
is steeper than that for L u , see IBerger et ail 119961 : 
Mendez et ail Il998bl : IMendez et all l2001h . Similarly, 
Mendez et all (|200lD use the same method to strengthen 



the identification of the single kHz QPO observed in 
the atoll source Aql X-l. To further investigate if this 
method could be used to identify the sharp kHz QPOs 
we report in Section 13 . 31 in Figure O we compare the 
energy depen dence of the kHz QPOs i n 4U 1608-52 
(IBerger et all 119961: IMendez et a!lll998bl: IMendez et all 
|200l|) and 4U 0614+09 ([Mendez et al.lll997h with that of 
Terzan 2. The data for Terzan 2 seem to fall in between 
those for Li and L u of 4U 1608-52 but shows a com- 
pletely different behavior than the data of 4U 0614+09. 
The fact that the rms amplitude of the upper kHz QPO 
in 4U 0614+09 and 4U 1608-52 are significantly different 
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(by up to a factor of 3) and that the data for Terzan 2 
fall in between those of Lt and L u in 4U 1608-52 show 
that the method does not lead to unambiguous results. 
Mean source luminosity, instantaneous luminosity and 
instantaneous QPO frequency may all affect QPO en- 
ergy dependence in addition to QPO type. 

5. SUMMARY 

1. We presented a detailed study of the time variabil- 
ity of the atoll source IE 1724-3045 (Terzan 2) 
which includes, for the first time, observations of 
this source in its island and banana states confirm- 
ing the atoll nature of this source. We find that 
the different states of Terzan 2 show timing be- 
havior similar to that seen in other NS-LMXBs. 
Our results for the extreme island st ate are consis- 
tent w ith t hose previously repo r ted inlBelloni et al.l 
(|2002fl and Ivan Straaten et alj (|2003fl . 

2. We report the discovery of kilohertz quasi-periodic 
oscillations (kHz QPOs). Although we do not de- 
tect two kHz QPOs simultaneously nor significant 
variability above 800 Hz, the detection of the lower 
and the upper kHz QPOs at different epochs and 
the power excess found at high frequencies (such as 
the case in intervals O or observation 90058-06-04- 
00) suggest that simultaneous twin kHz phenomena 
as well as significant variability up to ~ 1100 Hz 
(or more) is probable. 

3. By comparing the dependence of the rms ampli- 
tude with energy of kHz QPOs in the atoll sources 
4U 1608-52, 4U 0614+09 and Terzan 2, we show 
that this dependence appears to differ between 
sources and therefore cannot be used to unambigu- 
ously identify the kilohertz QPOs in either L u or 
Li, as previously thought. 

4. We studied the flux transitions or flares observed 
since February 2004 and from the source state 
changes observed we conclude that they are due 
to aperiodic changes in the accretion rate. 

5. State transitions between the extreme island state 
and the island state should be observable in obser- 
vations at the beginning or at the end of each flare. 
Given the relatively gradual and predictable tran- 
sitions, Terzan 2 becomes the best source known 
upto now to study such transitions. 
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rms 

(%) 


Component 
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130.8 ± 7.4 
8.5 ±0.20 
0.74 ±0.01 
0.512 ± 0.004 
0.131 ± 0.005 


0.21 ±0.07 

0.53 ±0.03 
6.5 ±0.9 


13.1 ±0.3 
18.5 ±0.1 
14.5 ±0.3 
3.9 ±0.2 
15.4 ±0.2 


-'-'low 

T h 

Llf 

Lb 
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169.2 ±8.6 
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0.92 ± 0.01 
0.750 ±0.006 
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0.31 ±0.07 
0.01 ±0.04 
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0.34 ±0.07 
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Lb 


1 nterva. 


H 


216.1 ± 2.2 
161.4 ± 17.0 
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1.21 ± 0.04 
0.93 ± 0.01 
0.179 ± 0.009 


17.9 ±8.4 
0.14 ±0.13 
0.10 ±0.07 
0.44 ± 0.04 
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3.1 ±0.5 
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* * * 
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T h 
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17.5 ± 2.4 
1.17 ±0.07 
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11 2 24 - 5 
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Lh 
Llf 
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Interval H 


131.7 ±15.8 

9.6 ±0.4 
0.892 ± 0.026 
0.634 ± 0.006 
0.133 ± 0.005 


0.3 ±0.1 

0.47 ± 0.04 
4.6 ±0.8 
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14.5 ±0.4 
4.5 ±0.4 
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1 


211.7 ±40.3 
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0.87 ±0.02 
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0.17 ±0.07 
0.47 ± 0.06 
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0.13 ±0.05 
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17.6 ±0.7 
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.1 
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* * * 
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68.5 ±7.2 1 1.8 ±1.1 1 5.5 ±0.9 1 *** 



Characteristic 
Frequency (Hz) 
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(%) 


Component 
ID 


Interval A 


156.9 ±6.6 
11.5 ± 0.2 
1.05 ±0.01 
0.820 ± 0.006 
0.171 ±0.005 


0.25 ± 0.06 
0.04 ± 0.03 
0.47 ± 0.02 
5.76 ± 0.87 
0.32 ± 0.06 


13.9 ±0.4 
17.6 ±0.3 

16.3 ±0.3 
3.3 ±0.3 

12.4 ±0.6 


'-'low 

L h 
Llf 
L h 


Interval H 


200.1 ± 11.6 
19.1 ± 1.5 
2.59 ±0.05 
0.40 ± 0.02 


0.55 ±0.09 
0.40 ± 0.04 


12.6 ±0.5 
13.2 ±0.2 
13.6 ±0.4 
12.2 ±0.2 


'-'low 

k 


Interval C 


148.4 ±6.8 
9.7 ±0.1 

0.86 ±0.01 
0.633 ±0.007 
0.158 ±0.006 
0.106 ±0.006 


0.27 ±0.05 

0.55 ± 0.02 
6.2 ± 1.2 

1.94 ± 1.02 


13.8 ±0.3 
18.6 ±0.1 
15.1 ±0.3 
3.1 ±0.3 
14.5 ±0.3 
3 4 +0.8 


'-'low 

L h 

J^LF 
Lb 

L v i 


Continued on next page 



(Umax = U centra l/FWHM), FRACTIONAL RMS (iN THE FULL PCA ENERGY 

band) and component identification (id) of the lorentzians fitted 

for Terzan 2. The quoted errors use Ax 2 = 1.0. Where only one 
error is quoted, it is the straight average between the positive and 

the negative error. " " : broad lorentzians whose quality factor 

was fixed to 0. 

**: This component is not 3ct significant. However, wi thou t this 

component the fit becomes unstable (see section 13.41 1 . 
* * *: For this component, no clear identification was possible. 
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